The inflationary prediction for primordial non-gaussianity 
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We extend the SN formalism so that it gives all of the stochastic properties of the primordial 
curvature perturbation £ if the initial field perturbations are gaussian. The calculation requires only 
the knowledge of some family of unperturbed universes. A formula is given for the normalisation 
/nl of the bispectrum of which is the main signal of non-gaussianity. Examples of the use of the 
formula are given, and its relation to cosmological perturbation theory is explained. 
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Introduction. The primordial curvature perturbation 
of the Universe, denoted here by £> is already present a 
few Hubble times before cosmological scales start to en- 
ter the horizon Its time- independent value at that 
stage seems to set the initial condition for the subse- 
quent evolution of all cosmological perturbations. As 
a result, observation probes the stochastic properties of 
C, which is found to be almost gaussian with an almost 
scale-invariant spectrum. 

According to present ideas £ is supposed to originate 
from the vacuum fluctuations during inflation of one or 
more light scalar fields, which on each scale are promoted 
to classical perturbations around the time of horizon exit. 
One takes inflation to be almost exponential (quasi de- 
Sitter spacetime) corresponding to a practically constant 
Hubble parameter and the effective masses of the 
fields to be much less than H*. This ensures that the 
fields are almost massless and live in almost unperturbed 
quasi de-Sitter spacetime, making their perturbations in- 
deed almost gaussian and scale-invariant. This automat- 
ically makes C almost scale-invariant, and can (though 
not automatically @>|j|) make it also almost gaussian. 

All of this is of intense interest at the present time, 
because observation over the next few years will rule 
out most existing scenarios for the generation of £, by 
detecting or bounding the scale-dependence and non- 
gaussianity of £. In this Letter we describe a general 
procedure for calculating the level of non-gaussianity, by 
means of the SN formalism 0, [f| . 

Defining the curvature perturbation. Perturbations of 
the observable Universe are defined with respect to 
an unperturbed reference universe, which is homoge- 
neous and isotropic. Its line element may be written 
as ds 2 = —dt 2 + a 2 {t)8ijdx' l dx^ defining the unperturbed 
scale factor a(t), time t, and the Cartesian spatial coor- 
dinates x. 

The curvature perturbation is only of interest after the 
universe has been smoothed on some scale (|) much 
bigger than the horizon To define it, one takes the 

fixed-i slices of spacetime to have uniform energy density, 
and the fixed- x worldlines to be comoving. The spatial 



metric is HHS0 

gij = a 2 (t)e 2 ^ x ) 7 ^,x) = a 2 (t,x) 7y (t,x) . (I) 

In this expression, 7y(<,x) has unit determinant, so that 
a volume of the Universe bounded by fixed spatial coor- 
dinates is proportional to the locally-defined scale factor 
a(i,x). In the inflationary scenario the factor 7^ just 
accounts for the tensor perturbation, but its form is ir- 
relevant here. According to this definition, £ is the per- 
turbation in In a. 

One can also consider a slicing whose metric has the 
form in Eq. without the £ factor, which we call the 
flat slicing. Starting from any initial flat slice at time i; n , 
let us define the amount of expansion N(t, x) = lnl ^4^- 
to a final slice of uniform energy density. Then |5j 



C(t, x) = 5N = N(t, x) - N (t) , 



(2) 



where No(t) = In 



a(ti, 



is the unperturbed amount of 



expansion. 

To make use of the above formalism we assume that 
in the super-horizon regime (aH 3> k), the evolution 
of the Universe at each position (the local evolution) is 
well-approximated by the evolution of some unperturbed 
universe 0, 0, @- This 'separate universe' assumption 
will presumably be correct on cosmological scales because 
these scales are so big 0. 

By virtue of the separate universe assumption, N (t, x) 
is the amount of expansion in some unperturbed uni- 
verse, allowing £ to be evaluated knowing the evolution 
of a family of such universes. For a given content of the 
Universe it can be checked using the gradient expansion 
El EiJ Ell method, but we do not wish to assume a 
specific content. 

The separate universe assumption leads also to local 
energy conservation, so that £ is conserved as long as the 
pressure is a unique function of the energy density. This 
consequence of the separate universe assumption was first 
recognised in full generality in Refs. 5, 10] (see also 
Ref. [a] for the case of inflation, Refs. lallJl for the case of 
linear perturbation theory, and Ref. [llj for a coordinate- 
free treatment). 
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The inflationary prediction. The evolution of the ob- 
servable Universe, smoothed on the shortest cosmological 
scale, is supposed to be determined by the values of one 
or more light scalar fields when that scale first emerges 
from the quantum regime a few Hubble times after hori- 
zon exit. Defined on a flat slicing, each field (pi at this 
epoch will be of the form </>j(x) = fa + 5<fii(x). 

Because quasi exponential inflation is assumed, and 
only light fields are considered, the perturbations 5<pi 
generated from the vacuum are almost gaussian, with 
an almost flat spectrum [l2| 



/ / 

2^ 



(3) 



Now we invoke the separate universe assumption, and 
choose the homogeneous quantities <j>i to correspond to 
the unperturbed universe. Then Eq. J5J for £ becomes 

C(t,x) - N(p{t), <Mx),^(x),---) -N(p(t), 0i,02, •••). 

(4) 

In this expression, the expansion N is evaluated in an 
unperturbed universe, from an epoch when the fields have 
assigned values to one when the energy density has an 
assigned value p. This expression 0, 0] allows one to 
propagate forward the stochastic properties of ( to the 
epoch when it becomes observable, given those of the 
initial field perturbations. 

Since the observed curvature perturbation is almost 
gaussian, it must be given to good accuracy by one or 
more of the linear terms (we use the notation N t i = 



and N , 



) 



(5) 



with the field perturbations being almost gaussian. In 
this Letter we include for the first time the quadratic 
terms 



(6) 



They may be entirely responsible for any observed non- 
gaussianity if the field perturbations are gaussian to suf- 
ficient accuracy. 

The bispectrum. The stochastic properties of the 
perturbations are specified through expectation values 
which, according to the inflationary paradigm, are taken 
with respect to the time-independent (Heisenberg pic- 
ture) quantum state of the Universe (to be precise, the 
quantum state of the universe before it somehow collapses 
to give the observed Universe). Focusing on £, we con- 
sider Fourier components, Ck = / d 3 k(,(t, x) exp(ik • x). 

The stochastic properties of a gaussian perturbation 
are specified entirely by the spectrum T>£, defined through 



(CkCk') = (2 7 r) 3 P c (fc) ( 5 3 (k 
From Eqs. © and © 



k') and V ( (k) ee ^P c {k). 



2^ 



(7) 



Non-gaussianity is defined through higher correlations. 
We consider only the three-point correlation. (The four- 
point correlation may give a competitive observational 
signature and can be calculated in a similar fashion 
|lj|.) It defines the bispectrum B^ through (CkCk'Ck") = 
(2ir) 3 B c (k, k', k")5 3 (k + k' + k"). Its normalisation is 
specified by a parameter /nl according to 0, 0] 
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/ NL (ft, k', fc")(P c (ft)P c (fc') + cyclic) 



(8) 



(In first-order cosmological perturbation the gauge- 
invariant gravitational potential $ during matter dom- 
ination before horizon entry is $ = — §Ci an d our 
definition of /nl coincides with the definition [l5| 
= 2/ NL (fc, ft', fc")[P # (fc)P$(fe') + cyclic]. At second- 
order these definitions of /nl differ \u\.) 

We shall take 7\ and /nl to be evaluated when cos- 
mological scales approach the horizon and ( becomes 
observable. Observation gives Vq = (5 x 10~ 5 ) 2 , and 
/nl| < E3- Absent a detection, this will eventually 
come down to roughly |/ NL | < 1 0. 

Ignoring any non-gaussianity of the 6<pi, our formula 
in Eq. JfjJ) makes /nl practically independent of the 
wavenumbers. Indeed, generalising the result found in 
Ref. we have calculated 

3 , E« NiNjNij Vc £ k N tij N Jk N M 

—/nl = — — —2 — ^ ln ( fci )^r 



In deriving this expression we used the spectrum (-^jr) of 
the field perturbations, and used Eq. JJJ to eliminate £/* 
in favour of Vq. As discussed in Ref. [ijj, the logarithm 
can be taken to be of order 1, because it involves the 
size ft -1 of a typical scale under consideration, relative 
to the size L of the region within which the stochastic 
properties are specified. Except for the logarithm, /nl is 
scale-independent if the field perturbations are gaussian. 
If only one d<f>i is relevant, Eq. © becomes 



C(t,x) =N i 64> i + -N ii (6<f> i ) 2 , 



(10) 



and because the first term dominates, Eq. © becomes 

(11) 



3 . _ 1 Nu 



In this case, /nl may equivalently be defined by writ- 
ing C = Cg - f/NLCg, where C g is gaussian. 

To include the possible non-gaussianity of the 5<pi, we 
define the bispectra B^ of the dimensionless field per- 
turbations (2ir / H*)8<pi and their normalisation fijk, in 
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exactly the same way that we denned Bq and /nl- These 
bispectra add the following contribution to /nl in Eq. © 



A/i 



NL — ^75 V c 



(12) 



The /ijfe , generated directly from the vacuum fluctuation, 
will depend strongly on the wavenumbers. 

Cosmological perturbation theory. In the super-horizon 
regime the non-linear theory || that we have used is a 
complete description. The basic expression in Eq. (0} 
is non-perturbative, giving £(i, x) in terms of the initial 
fields and the expansion of a family of unperturbed uni- 
verses. The second-order expansion in Eq. © is a matter 
of convenience. As we shall see it seems to be adequate 
in practice, but Eq. would still be applicable if the 
expansion converged slowly or not at all. 

Cosmological perturbation theory (CPT) is completely 
different. It is applicable both inside and outside the hori- 
zon, being at each instant a power series in the perturba- 
tions of the metric and the stress-energy tensor, together 
with whatever variables are needed to completely specify 
the latter and close the system of equations. During in- 
flation these variables are the components of the inflaton, 
while afterwards they may involve oscillating fields and 
a description of the particle content. First-order CPT 
is usually adequate and can describe non-gaussianity at 
the level /nl 3> 1, which has to be generated by the 
second-order term in Eq. Second order CPT is gen- 
erally needed only to handle non-gaussianity at the level 
|/nl| - 1- 

Quantised CPT is needed to calculate the stochastic 
properties of the initial field perturbations S(f>i , which are 
the input for our calculation. The slow-roll spectrum in 
Eq. comes from the first-order calculation. The bis- 
pectrum is a second-order effect and has, in the context 
of slow-roll inflation, been calculated in Refs. 0, H^| . 
It is shown elsewhere 01 that |A/nl| C 1 in this 
case. Higher correlators have not been calculated yet 
and would give an additional contribution to Eq. © 
which presumably is also negligible. Exotic non slow-roll 
models [2l] can make |A/nl| ^ I, but from now on we 
set A/nl = 0. 

In the regime aH 3> k, perturbation theory must be 
compatible with Eq. In particular, the non-local 

terms, present at second order for a generic perturbation, 
must be absent for ( (see also Ref. [§2]). Finally, CPT 
is needed to evolve the perturbations after horizon entry, 
but that is not our concern here. In the following, we 
apply our formalism to calculate /nl in various cases and 
compare it with the CPT result where that is known. 

A two- component inflation model. As a first use of 
Eq. © we consider the two-component inflation model of 
Kadota and Stewart [13] , estimating for the first time the 
non-gaussianity which it predicts. The model works with 
a complex field which is supposed to be a modulus 



with a point of enhanced symmetry at the origin. Writ- 
ing $ = |$|e , the tree-level potential has a maximum 
at $ = and depends on both |$| and 9. A one-loop 
correction turns the maximum into a crater and inflation 
occurs while $ is rolling away from the rim of the crater. 
The curvature perturbation is supposed to be constant 
after the end of slow-roll inflation. For 9 <C 9 Cl with 9 C 
being a parameter of the model, it is found that N oc | ^ | . 

Through the first term of Eq. /nl — ■§- which is too 
small ever to be observed. 

The curvaton model. In the curvaton model [24| 
(see also Ref. 25]) the curvature perturbation ( grows, 
from a negligible value in an initially radiation domi- 
nated epoch, due to the oscillations of a light field a 
(the curvaton) around the minimum of its quadratic po- 
tential V a (t,x) = imJ(T 2 (t,x), where m a is the cur- 
vaton effective mass. Due to the oscillations, the ini- 
tially negligible curvaton energy density redshifts as 
p CT (i,x) w |mj^(t,x) cx a~ 3 (£,x), where a a represents 
the amplitude of the oscillations. Meanwhile the radia- 
tion energy density p r redshifts as a~ 4 . Soon after the 
curvaton decay, the standard Hot Big-Bang is recovered 
and £ is assumed to be conserved until horizon reentry. 

To calculate /nl using Eq. © we first realise that <r» 
(the value of a a few Hubble times after horizon exit) 
is the only relevant quantity since the curvature pertur- 
bation produced by the inflaton, and imprinted in the 
radiation fluid during the reheating process, is supposed 
to be negligible. Thus, Eq. (f 1 1 1> applies. Second, we can 
redefine N as the number of e-folds from the beginning 
of the sinusoidal oscillations to the curvaton decay. This 
is because the number of e-folds from the end of inflation 
to the beginning of the oscillations is completely unper- 
turbed as the radiation energy density dominates during 
that time. Thus, N is now a function of three variables 



N(p 

dec j Po 



1 



hi 









3 



P<T d e 



(13) 

where g = a osc is the amplitude at the beginning of 
the sinusoidal oscillations, as a function of <7„. Here 
the curvaton energy density just before the curvaton de- 
cay Pa dcc is expressed in terms of the total energy den- 
sity pdec at that time, the total energy density at the 
beginning of the sinusoidal oscillations p sc, and g by 

[g{cr,)} 2 ( Pde p3r deC ) 3/4 - Mter evaluatin S 
sc , we obtain 



1 2 r 
2 m o 



Po-dec ~ 



■&k = 9' 3£> at nxed Pd, 



and p . 



N, 



2 g' 
-r- 



(14) 



3 9 

being p TdFa the radiation energy 



where rs. , 
density just before the curvaton decay, so that 



2tt 



-JV., 



H*r g' 
3tt g ' 



(15) 
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in agreement with first-order cosmological perturbation 
theory Q. Differentiating again we find from Eq. (|f f I) 



NL 



6 N* 



5 5 
- + -r 

3 6 



^(1 + ^) 'I* 



4r 



which nicely agrees with the already calculated /nl usm 
first- and second-order perturbation theory (see Refs 

ilH). 

Another two- component model. Finally we consider 
the two-component inflation model of Refs. 0, H3|- 
For at least some number N of e-folds after cos- 
mological scales leave the horizon, the potential is 

V = V (l + 3»?*;£r + 5»7cr;*r), with the first term 
dominating, r\^ and r\ a being the slow-roll 77 parameters, 
and mp being the reduced Planck mass. The idea is to 
use Eq. © to calculate the non-gaussianity after the iV 
e-folds which, barring cancellations, will place a lower 
limit on the observed non-gaussianity. 

The slow-roll equations give the field values <j>(N) and 
o~(N), in terms of those obtaining just after horizon exit; 
(j)(N) = 0exp(-iV?70) and a(N) = a exp(-Nrj a ). This 
gives V(N, 4>, o) and allows us to calculate the derivatives 
of N with respect to <f> and a at fixed V. Focusing on the 
case a = considered in Ref. |27|, we find 



_ 84> 7)4, / 8<t> \ | V* e 2 N (r,+ -r,„) f _^_\ ' 



(17) 



in agreement with the second-order perturbation calcu- 
lation of Ref. [2^1 . If the observed £ has a non-gaussian 
part (a equal to the last term of Eq. fTTjl and a gaus- 
sian part generated mostly after inflation, one can obtain 
|/nl| > 1 by choosing ^ > 0.26, n a = ^,N = 70, and 

This model was studied originally |3j, |27j using 
a second-order perturbation expression for the time- 
derivative of eHmp^ir, with e being the e slow- roll pa- 
rameter. This expression disagrees with ours |29| through 
the appearance of non-local terms, though the order of 
magnitude is similar [sol ]. 
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